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Abstract 


This paper presents a comprehensive procedure to 
predict both the vibration and noise generated by a gear 
transmission system under normal operating conditions. 
The gearbox vibrations were obtained from both numeri- 
cal simulation and experimental studies using a gear 
noise test rig. In addition, the noise generated by the 
gearbox vibrations was recorded during the experimental 
testing. A numerical method was used to develop linear 
relationships between the gearbox vibration and the gen- 
erated noise. The hypercoherence function is introduced 
to correlate the nonlinear relationship between the fun- 
damental noise frequency and its harmonics. A numeri- 
cal procedure was developed using both the linear and 
nonlinear relationships generated from the experimental 
data to predict noise resulting from the gearbox vibra- 
tions. The application of this methodology is demon- 
strated by comparing the numerical and experimental 
results from the gear noise test rig. 


Nomenclature 

A rotor modal displacement of (X, 0 X ) 

A c casing modal displacement of (X c , 0^) 

B rotor modal displacement of (Y, 0 y ) 

B c rotor modal displacement of (Y c , 0 cy ) 

[C b ] bearing damping matrix 

[C ] casing structure damping matrix 

[Cj i*] T [CJ [*] - 

[cj [*cl T ICJ w 

[c c ] [* e ] [cj (*J 

D rotor modal displacement of Z 

D t rotor modal displacement of 0 t 

D c casing modal displacement of (Z c , 0 CE ) 


(F(t)} 

(F G (t)} 

(F c (t)} 

{F.(t)} 

[g a ] 

[G v ] 

(g a 1 

[Gy] 

[K b ] 

[K c ] 

[KJ 

[k a 1 

[KJ. 

[K a ] 

[K b ] 

K b x>K b y 

[M] 

! M cl 

w 

W c 

x 

x c 

y 

Vc 




external and mass-imbalance excitations 
gear force 

force acting on casing structure 
shaft bow force 
rotor angular acceleration 
gyroscopic 

W T I g aI 1*1 
!*1 t [g.1 1*1 

bearing stiffness 
casing structure stiffness 
shaft bow stiffness 
l/2{[K b J + [K by ]} 

i*] T iK b i m 
[*] T [K a ] {§] 

\*f [KJ l* c ] 

bearing stiffness in x and y directions 

mass matrix of rotor 

mass matrix of casing structure 

generalized displacement vector of rotor 

generalized displacement vector of casing 

rotor x-displacement 

casing x-displacement 

rotor y-displacement 

casing y-displacement 

rotor z-displacement 

casing z-displacement 

friction coefficient between the gear teeth 
surface 
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$ 

*x 

$ x<> 


v x* v y iV t 

*cx.*cy> 


w 

w. 


angle of orientation 
rotor orthonormal mode shape 
casing orthonormal mode shape 
rotor tran slat ion al mode shape 
rotor rotational mode shape 
rotor rotational motion 
casing rotational motion 

critical speed of rotor 
critical speed of casing 


Introduction 


The ever increasing power- to- weight ratio require- 
ments of newer transmission systems often result in 
greater noise and vibration at both the gear mesh and 
the gearbox structure. Large vibrations in gear transmis- 
sion systems result in higher gear tooth wear rates and 
possible tooth root crack formation which leads to pre- 
mature gear failure. In addition, excessive noise pro- 
duced by a gear transmission in an aircraft causes crew 
fatigue, strained communication, and possible hearing 
damage for crew and passengers. In order to assure a 
quiet, smooth, and safe operation of a gear transmission 
system it is necessary to understand the dynamics of the 
system and the noise generated by the system under var- 
ious operating conditions. 

A large amount of work is reported in the literature 
on analyzing the localized effects of gear tooth interac- 
tions. In the area of experimental investigations of 
localized gear tooth interactions, a significant amount of 
work was performed over the last 10 years (Akin, 1989; 
Choy, 1989; Oswald, 1987; Krantz, 1992; Lewicki, 1992; 
El-Bay oumy, 1989; and Rebbechi, 1991). There is much 
less information available in the literature on the study 
of the global dynamics of a transmission system. There 
have been some studies on the vibration analysis of a 
single-stage gear (August, 1982; Choy, 1988; Mark, 
1982; and Boyd, 1989) and on multistage gear dynamic 
analysis (Choy, 1991, 1993; David, 1987; and Ozguven, 
1988). The experimental study of gear system vibrations 
is somewhat limited (Choy, 1991; Lim, 1989; and 
Oswald, 1992). Reported research is even more limited 
in the areas of noise analysis and noise prediction in gear 
transmission systems. Some work on predicting and veri- 
fying the noise field of a gear transmission system has 
only recently been reported (Seybert, 1991 and Oswald, 
1992). 


In response to the importance of the application and 
the limited amount of work performed in the area of 
gearbox noise and vibration research, a program was 
started with the main objective of developing the meth- 
odology required to predict gear transmission system 
noise and vibration. This paper describes the global dy- 
namic model that was developed to simulate the noise 
and vibration of a multistage gear system. Results pre- 
dicted by the numerical model are compared to experi- 
mental results obtained from the gear noise test rig at 
NASA Lewis Research Center. 


Analytical Procedure for Gear System Vibration 

To achieve the numerical solution, the generalized 
equations of motion developed in a previous paper 
(Choy, 1993) presents the numerical procedure using the 
modal synthesis method. Using the undamped modes for 
both the rotor-bearing-gear system and the gearbox 
structure, the generalized equations of motion for the 
system are transformed into the corresponding modal co- 
ordinates. This procedure significantly reduces the 
degrees of freedom of the system, and provides an esti- 
mate of the modal excitations at the various modes. 


The modal transformation procedure can be devel- 
oped by assuming that the generalized displacements can 
be represented by a linear combination of the modal 
characteristics, (Choy, 1991, 1993) as 


{W} 


f* J {A} 
fcrfKA) 
[*yl (B> 

MW 
f»J fl>> 

[*tl (DJ 


(1) 


where A, B, D, and D t are the time-dependent modal 
excitation functions, and [$ x ], [$ y ], [$ y 0]> f®*]i 

and are the matrices containing the orthonormal 
modes of the corresponding rotor system. Similarly, for 
the gearbox housing structure, 
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[$cy]{B c } 
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(2) 


using 


{Z} 


{A} 

{B} 

K 

{D} 

(DJ 


{ZJ 


{A c } 

{B c } 

{D c } 


(3) 


The modal equations for the rotor system (Choy, 1991, 
1993) can be written in modal coordinates as: 

[i] {z} = [n v ]{z} - [G A ]{Z} - [C b ]{z} 

+ (K b -K A l{Z}-[$] T [C b ][$ c ]{Z c } 

- |w 1 2 ]{Z} + [*1 T [K b ][* c ]{Z c } 

+ [$] T {F(t) + F G (t)+F g (t)} 

(4) 

and for the casing 


[I] {ZJ = [CJ {ZJ - [ W J] {ZJ - [K b ] {ZJ - [C b ] {ZJ 

+ [*J T (K b ] 1*1 {Z} + [# C ] T [c b ] [$] {Z> + [$ C ] T {F c (t)} 

(5) 

Using the appropriate initial conditions for the modal 
velocity {Z} and displacement {Z}, the modal accelera- 
tions {Z} can be solved. The solution procedure can be 
summarized as follows: 


2. Use variable time-stepping integration scheme to 
evaluate modal velocity and displacement at next time 
step. 

3. Calculate velocity and displacement in generalized 
coordinates. 

4. Calculate nonlinear gear mesh forces due to rela- 
tive motion between the two gears. 

5. Calculate bearing forces due to relative motion 
between the rotor and the casing. 

6. Calculate modal bearing and gear forces. 

7. Repeat from step number 2. 


Experimental Studies 

Test results from the gear noise test rig (Fig. 1(a)) at 
NASA Lewis were used in this study to refine and vali- 
date the global dynamic model developed. Three phases 
of experimental studies were performed: (l) the determi- 
nation of static modal characteristics using the impulse- 
response technique; (2) the transient vibration of both 
rotor and the gearbox structure during operation; and 
(3) measurement of the noise generation at the gearbox 
surfaces. An outline of the experimental studies is pre- 
sented below: 


Evaluation of Experimental Modal Characteristics 

The gear transmission system used in this test is 
shown in Fig. 1(b). The gear system consists of a pair of 
identical parallel axis spur gears supported by rolling 
element bearings. The top and four sides of the gearbox 
were divided into a grid of 116 node points. The transfer 
functions of the system were found from the impulse- 
response technique by using an impact hammer. Three 
accelerometers, one on the surface in the x- direct ion, 
another in the y-direction, and the last in the 
z-direction, were used to obtain the responses from the 
impacts. Natural frequencies and their corresponding 
mode shapes excited during the experiment were evalu- 
ated using a two-channel dynamic signal analyzer with 
PC-based modal software. The modal characteristics 
found during the experiment were compared with those 
generated by the numerical model for verification 
purposes. 


Experimental Investigations of System Dynamics 


1. Input initial conditions to calculate modal 

accelerations. 


The test rig is powered by a 150-kW (200-hp) vari- 
able speed electric motor. An eddy-current dynamometer 
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is used to absorb power at the opposite end. During this 
test, accelerometers were placed on the top surface of the 
gearbox to monitor the vibrations of the gearbox struc- 
ture. In addition, two sets of noncontacting eddy-current 
proximitors monitored motion of the gear shafting. The 
vibration signals of the gearbox surface and the shafts 
were recorded over a range of operating speeds from 
3000 to 5000 rpm. A two-channel dynamic signal anal- 
yzer computed the vibration frequency spectra which are 
displayed in a waterfall diagram. The lateral vibration 
frequency spectra of the gears were also computed. The 
results of this phase of the experimental work were com- 
pared with those from the model for verification of the 
numerical procedure and improvement of the modeling 
of the gear transmission system. 

Investigations of Gearbox Noise 

Two acoustic microphones were used in this experi- 
ment. Both microphones were placed over the top cover 
of the gearbox with one in the vertical direction and the 
other at a 45° angle from the vertical. While rotor and 
gearbox vibration are recorded as described in the last 
section, noise was obtained by the microphone system 
at various running speeds. A correlation between the 
gearbox vibrations and the noise level was performed in 
the frequency domain using linear transfer and cross- 
correlation functions. The nonlinear hypercoherence 
function (Jong, 1986) was also used to generate a rela- 
tionship between the noise and the vibration data in 
order to develop a noise prediction algorithm for the 
global dynamic model. 

Correlation and Benchmarking of System Dynamics By 

Numerical Procedures with Experimental Results 

The modal hammer excited many vibration modes of 
the gearbox structure. Eight of these modes in the fre- 
quency range 0 to 3000 Hz seem to dominate the trans- 
fer function spectra. These eight modes, shown in Fig. 2, 
were chosen to represent the vibration of the system. In 
order to produce an accurate analytical simulation of the 
test gearbox, a similar set of modes was generated using 
a NASTRAN finite element model of the gearbox struc- 
ture. Out of 25 modes in the analytical model in the 
0 to 3000 Hz frequency range, the eight dominant modes 
(similar to those from the experimental studies) were 
used to represent the gearbox dynamic characteristics. 
These analytically simulated modes are shown in Fig. 3. 
Table I lists the natural frequencies of the simulated 
modes. The calculated natural frequencies are within 
5 percent of the measured modes. The three-dimensional 
analytical mode shapes (Fig. 3) are very similar to the 
experimental mode shapes (Fig. 2). The correlation of 
the results between the analytical model and the experi- 


mental measurements confirms the accuracy of the dy- 
namic representation of the test gearbox using only a 
limited number of modes. In addition, this correlation of 
the experimental modal characteristics with the numeri- 
cal model helps to reduce the number of modes being 
used in the numerical simulation to achieve an accurate 
result. 

A slow roll of the rotor-gear assembly revealed that 
a substantial residual bow was present in the rotor as 
shown by its large orbital motion given in Fig. 4. Fig- 
ure 4(a) represents the orbit of the driver rotor at the 
gear location and Fig. 4(b) represents the orbit of the 
driven rotor. The circular orbit in the driver rotor at 
slow roll represents the residual bow deformation of the 
rotor. The elliptical orbit in the driven rotor is due to a 
combination of the residual bow effects and the vertical 
gear force from the torque of the driving gear. In order 
to analytically simulate the influence of this effect, a 
residual bow of 0.05 mm (2 mil) was incorporated into 
the numerical model. Figure 5 represents the frequency 
spectra of the driven gear vibration generated from both 
the experimental and the numerical procedures. The ex- 
cellent agreement in both results further verifies the 
validity of the numerical study. 

In order to predict the noise of the transmission 
system during operation, the vibration of the gearbox 
must be accurately modeled. Figures 6 and 7 present the 
time domain vibration signal and corresponding fre- 
quency spectra of the signal for the gearbox top surface 
from both experimental and numerical studies. Fig- 
ure 6(b) shows the analytically predicted vibration sig- 
nal (acceleration) at the top of the casing structure at an 
operating speed of 3000 rpm while the corresponding ex- 
perimental results are given in Fig. 6(a). Note that the 
amplitudes and general shapes of the two time signals 
are very similar. The comparison of their frequency 
contents are given by the frequency spectra in Fig. 7. 
Figure 7(a) shows the frequency components of the 
measured vibration at the top of the gearbox. Note that 
the major vibration components occur at the mesh fre- 
quency of 1250 Hz and at 3-tunes (3750 Hz) and 4-tunes 
(5000 Hz) the mesh frequency. The very large amplitude 
in the frequency components at 3750 and 5000 Hz are 
due to the excitation of the gearbox natural frequencies 
near the 3-times and 4-times multiples of the mesh fre- 
quency. The fundamental and the 2-times mesh fre- 
quency component are substantially smaller due to the 
lack of any gearbox natural frequencies near 1250 and 
2500 Hz. In addition, frequency components of substan- 
tial magnitude are also noticed at 700 and 2900 Hz, 
which are also close to the natural frequencies of the 
gearbox. Figure 7(b) depicts the results from the numer- 
ical simulations of the system. The numerical results 
correspond closely to those from the experimental study 
with the exception of the 4-times mesh frequency 
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component. The magnitude and the general distribution 
of the frequency content of the spectrum is very similar 
to the experimental results. 

Correlation of Noise and Vibration 

In order to predict the noise generated in a gear 
transmission system, a two-phase procedure was devel- 
oped to correlate the noise and vibration data obtained 
during experimental studies: (l) linear correlation 

through transfer functions, and (2) nonlinear correlation 
through hypercoherence functions. The relationships 
derived in this study were used to predict gearbox noise 
from the numerical simulation of the gearbox vibrations. 
A description of the procedures used is presented in the 
following paragraphs: 

Correlation of Experimental Noise and Vibration Data 

In this experiment, a new set of two identical 25- 
tooth gears were installed in the gearbox shown in 
Fig. 1(b). The vibration of the gearbox surfaces was 
monitored by accelerometers, and the noise data was 
obtained by microphones, as described in the previous 
section. The frequency and power spectra of both the 
vibration and noise signals at 3000 rpm running speed 
are shown in Fig. 8. Note that in the vibration spec- 
trum, a small frequency component exists at the mesh 
frequency of 1250 Hz with two very large components at 
3-times the mesh frequency (3750 Hz), and 4-times mesh 
frequency (5000 Hz). In the noise spectrum, while fre- 
quency components of considerable magnitudes are 
observed at the mesh frequency of 1250 Hz, and at 2 and 
3-times the mesh frequency (2500 and 3750 Hz), a large 
amplitude frequency component is also noticed at 
400 Hz. This 400-Hz frequency component does not coin- 
cide with any of the major vibration frequencies at the 
top surface of the gearbox. A closer examination of the 
system dynamic characteristics reveals that the 400-Hz 
frequency is an excitation of one of the natural frequen- 
cies of the rotor system, which has a very small ampli- 
tude in the gearbox vibration spectrum. A logical 
explanation is that noise is being produced at the rotor- 
gear system at 400 Hz and is not being transferred to 
the vibration of the gearbox surface. This phenomenon 
is further verified by the very large amplitude of the 
vibration-to-noise transfer function at 400 Hz as given 
in Fig. 9(a). A linear coherence function between the 
noise and vibration signals is also given in Fig. 9(b). 

From the experimental data, it can be observed that 
in the frequency domain, most of the noise components, 
other than the 400-Hz contributions, exist at multiples 
of the mesh frequency. The nonlinear hypercoherence 
function was used for the second phase of this study. 


The hypercoherence function (Jong, 1986) establishes a 
nonlinear relationship between the basic frequency com- 
ponent and its multiple harmonics. Figure 10 shows the 
hypercoherence functions of the noise signals for the 
3000-rpm case. The magnitudes of the harmonic compo- 
nents form the nonlinear relationship of each component 
to the fundamental gear mesh frequency. These relation- 
ships provide a picture of the relative contribution of the 
higher harmonics with respect to the fundamental fre- 
quency. These hypercoherence functions were used on 
the numerical simulation data to examine the accuracy 
of the noise prediction. 

In order to predict the noise level of the transmission 
system during operation, the vibrations of the gearbox 
structure were simulated using the modal synthesis 
method. The transfer function for microphone number 1 
(vertical to the plate), as given in Fig. 9(a), was used to 
predict the noise spectrum of the system (Fig. 11(b)). 
Figure 11 shows the comparison of the results of the 
experimental noise (Fig. 11(a)) with the predicted noise 
(Fig. 11(b)) for the vertical microphone. Although all 
the major frequency components excited in both cases 
are very similar, one can see that the predicted noise 
spectrum (Fig. 11(b)) has a somewhat higher amplitude 
at mesh frequency (lx) and. a very small 4-times mesh 
frequency (4x) component as in the vibration spectrum. 

Another way of examining the validity of the pre- 
dicted noise signal is through the nonlinear relationships 
derived by the hypercoherence function. Figure 12 shows 
the comparison of the hypercoherence functions from 
both experimental and numerical solutions for micro- 
phone number 1. Note that both the experimental and 
numerical results show a very small 2-times component, 
with a substantial amplitude at 3, 4, and 5 times funda- 
mental frequency. However, a 50-percent coherence is 
found in the 3, 4, and 5 times components from numeri- 
cal predictions (Fig. 12(b)) while only about a 25- 
percent coherence is seen from the experimental results 
(Fig. 12(a)). This difference is possibly due to the 
limitations of the numerical model in simulating various 
details in the gearbox transmission system. An explana- 
tion of the much higher (nearly double) coherence ampli- 
tudes in the numerical study is that only the major 
frequency-related vibration factors are included in the 
model, which results in a predicted vibration signal con- 
sisting mostly of frequency components that are closely 
related to the fundamental frequency. 

Conclusions and Summary 

A methodology was developed using the modal syn- 
thesis technique to simulate the global dynamics of a 
gear transmission system. Good correlation has been 
obtained using results from the global dynamic model 
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and measured dynamic characteristics from the NASA 
gear noise test rig. A numerical relationship was also 
developed from the experimental data to correlate noise 
and vibrations in the frequency domain. This relation- 
ship was applied to the vibration simulations to predict 
the noise generated by the gearbox. Results of the pre- 
dicted noise were compared to the experimentally meas- 
ured data linearly using the frequency spectra and 
nonlinearly using the hypercoherence function. Some 
conclusions of this study can be summarized as follows: 

1. The dynamics of a gear transmission system 
coupled with the vibrations of the gearbox structure 
during operation can be simulated within 5 percent 
of the modal frequencies using the modal synthesis 
technique. 

2. Predicted vibrations at the housing surface using 
the global dynamic model are in good agreement with 
measured values. 

3. Both the frequencies and magnitudes of the pre- 
dicted noise components are similar to those obtained 
from experimental study. 

4. The hypercoherence function was found to be a 
valuable tool in comparing the nonlinear relationships in 
the predicted noise to the corresponding nonlinear rela- 
tionships in the measured noise. 
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TABLE I.— COMPARISON OF EXPERIMENTAL 
MEASURED AND ANALYTICAL MODELED 
NATURAL FREQUENCIES 


Experimental, 

Hz 

Analytical, 

Hz 

Difference, 

percent 

658 

658 

0 

1049 

1006 

-4.1 

1710 

1762 

3.0 

2000 

2051 

2.6 

2276 

2336 

2.6 

2536 

1536 

0 

2722 

2752 

1.1 

2962 

3012 

1.7 
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Dynamometer-^ 



(a) Schematic of NASA gear noise rig. 



(b) Test gearbox. 

Figure 1 .—Gear noise test rig at NASA LeRC. 
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Mode: 1 
Frequency, 658.37 Hz 


Mode: 2 
Frequency, 1 048.56 Hz 


Mode: 3 
Frequency, 1 709.77 Hz 


Mode: 4 


Frequency, 1999.95 Hz 


£ 
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Mode: 5 
Frequency, 2275.69 Hz 


Mode: 6 


Frequency, 2535.77 Hz 




Mode: 7 Mode: 8 

Frequency. 2722.1 6 Hz Frequency, 2961 .71 Hz 

Figure 2. — Experimentally determined mode shapes of the gear box (0 to 3000 Hz) 
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Mode: 1 

Frequency = 658 


Mode: 2 

Frequency = 1 006 



Mode: 5 Mode: 6 

Frequency = 2336 Frequency = 2536 



Mode: 7 Mode: 8 

Frequency = 2752 Frequency = 301 2 

Figure 3. — Analytically determined mode shapes of the gear box (0 to 3000 Hz). 
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Shaft #1 



Shaft #2 



Displacement, X-direction, In. 

(a) Driving rotor orbit. (b) Driven rotor orbit. 

Figure 4.— Orbital motion of rotor during slow roll. 
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Frequency, Hz 

X-direction vibration Y-direetion vibration 

(b) Analytical. 

Figure 5.— Frequency spectra of the driving rotor vibration. 
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Noise amplitude Vibration amplitude -g. Amplitude, volts Amplitude, g 



(a) Frequency spectrum for vibration 



phone number 1. 


2000x103 



(c) Power spectrum for vibration. 



(d) Power spectrum for noise from microphone 
number 1 . 

Figure 8.— Frequency and power spectrum at 
the operating speed of 8000 rpm. 




0 2000 4000 6000 8000 10 000 


Frequency, Hz 
(b) Coherence function. 

Figure 9. — Functions between noise (microphone number 1) and 
vibration at the operating speed of 3000 rpm. 
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Harmonics of mesh frequency 

Figure 10. — Hypercoherence function for the noise signal from 
microphone number 2 at 3000 rpm. 
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(b) Predicted noise. 

Figure 11 . — Frequency spectra results from microphone number 1 
at the operating speed of 3000 rpm. 


Frequency harmonic 
(b) Numerical predictions. 

Figure 1 2. — Hypercoherence function of the noise signal from 
microphone number 1 at 3000 rpm. 
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